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Using a central limit theorem for arrays of interacting quantum systems, we give analytical ex- 
pressions for the density of states and the partition function at finite temperature of such a system, 
which are valid in the limit of infinite number of subsystems. Even for only small numbers of 
subsystems we find good accordance with some known, exact results. 

Keywords: quantum statistical mechanics, spectral densities, partition sums, quantum central limit theorem 









X3 

o 

> 

O 

o 

O 

o 

-(— > 



I 

o 
o 



X 



I. INTRODUCTION 

In order to study a quantum system in full detail, its 
Hamiltonian needs to be diagonalized. With increasing 
dimension of the Hilbert space, the diagonalization of an 
operator becomes a very tedious task. For lattices or ar- 
rays of interacting subsystems or particles, the dimension 
of the Hilbert space scales as m" , where m is the dimen- 
sion of the Hilbert space of one particle and n is the 
number of particles. Thus, apart from some exceptions 
|3J , it is even numerically impossible to exactly determine 
the eigenvalues of those models. 

Fortunately, considerable understanding can be ob- 
tained already from functions of the eigenvalues with- 
out knowing each individual eigenvalue. For example all 
thermodynamical quantities of a system are determined 
by its partition function [a, |3, [l3 • 

In this paper, we present a novel approach, that allows 
to derive analytical expressions for all quantities whose 
operator is a function of the Hamiltonian for chains or 
lattices of very many interacting particles [ll l^, U^ ■ 

Our approach is based on a central limit theorem 
|l3lll4l | for quantum systems with nearest neighbor inter- 
actions. It is also valid if each subsystem does not only 
interact with its nearest neighbors but with a fixed, finite 
number of neighbors. 

Systems of that interaction topology play a central role 
in condensed matter theory 0,113 • Most models, which 
are currently used to describe strongly interacting elec- 
trons, belong to this category. 

In the analysis of the thermodynamics of systems of 
interacting particles, potential phase transitions are of 
central interest. Being based on a central limit theorem, 
our approach becomes more exact with increasing num- 
ber of subsystems. It is precisely this limit of infinite 
number of subsystems, that is relevant for the study of 
the possible phase transitions [l^, [lOI . 

The paper is organized as follows: In section ^ we 



present the class of models we address an introduce the 
notations we use. Section IIIII contains the central limit 
theorem, our approach is based on. The proof is not 
displayed here as it already appeared in a previous pub- 
lication [13. In the following two sections . II VI and Ivl we 
give the analytical expressions for the density of states 
and the partition function, which are a straight forward 
application of the theorem of section IIIII In the next 
section IVTl we numerically evaluate these expressions for 
an Ising spin chain and compare our results with the re- 
sults of an exact diagonalization. Section IVlII contains a 
discussion of some limits and problems of our approach. 
Finally section rVIIII summarizes the results. 



II. MODEL AND NOTATION 

We consider a chain of quantum systems with next 
neighbor interactions. The entire system is described by 
a Hamiltonian H which is a linear, self-adjoint operator 
on a separable, complex Hilbert space J^. The Hilbert 
space J^ is a direct product of the Hilbert spaces of the 
subsystems. 



^-n 



y^^fj,: 



and the Hamiltonian may be written in the form. 



H 






(1) 



(2) 



with 



lp.,ti+l'' 



'Electronic address; Imichael . hartmann@dlr . del 



j-[]-«in-(M+l) 

(3) 

where i/^ is the proper Hamiltonian of subsystem /i, and 
Ifj..fj.+i the interaction of subsystem /i with subsystem 
^ + 1. I is the identity operator. We chose periodic 
boundary conditions /„.„+i = In,i- 

Let Eu; be the eigcnenergies and, using the Dirac no- 
tation Q, let {\ip)} be an orthonormal basis of ,^ con- 
sisting of eigenstates of the total system. 

H\ip) = E^\ip) with {ip\ip') = S^^, , (4) 



where b^p^pl is the Kronccker delta. We denote by |a) the 
product states 



n 

\a)= n®i«M>> 



(5) 



built up from some eigenstates \a^^ of each subsystem 
Hamiltonian iJ^j, iJ^|a^) = E^\a^. Let Ea denote the 
sum of all _E^, Ea = X]u=i^m- ^^^ states \a) are as- 
sumed to form an orthonormal basis of J^if. We further- 
more define, 



Ea = {a\H\a), 

Al = (alH^la) ~ (a\H\a)\ 



and the operator 



H-Ea 



Z has the same eigenstates as H, 



(6) 

(7) 



(8) 



(9) 



where the z^ denote its eigenvalues. Note that H and 
therefore Z, E^p, z^p, Ea and Aq as well as the two bases 
{|(^)} and {\a)} depend on n. 

The measure of the quantum mechanical distribution 
of the eigenvalues of Z in the state \a), is given by the 
usual formula, 

F,(z^g[zi,Z2])= ^ |(a|^)p, (10) 

{\'p):zi<z^<Z2} 

where the sum extends over all states \tp) with eigenvalues 
in the respective interval. 



III. CENTRAL LIMIT THEOREM 

If the operator H and a state \a) satisfy 

Al>nC (11) 

for all n and some C > and if each operator 7i^ is 
bounded, i.e. 



(xI^mIx) < C 



(12) 



for all normalized states |x) G =^ and some constant C", 
then the quantum mechanical distribution of the eigen- 
values of Z in the state \a) converges weakly to a Gaus- 
sian normal distribution: 



lim Faiz^ G [21,22]) = 

n — ^00 

for all —00 < zi < Z2 < 00. 



exp (- z2/2) 



/2tt 



(13) 



The rigorous proof of this theorem is given in [15J. 
Note that the theorem also holds for lattices of inter- 
acting quantum systems in arbitrary dimension. Fur- 
thermore the interaction need not be limited to nearest 
neighbors. As long as each particle only interacts with a 
fixed number of neighbors, the theorem holds. 

If each subsystem H^ has an infinite energy spectrum, 
condition ifT^ is not fulfilled. The constant C in ifT^ 
however may be chosen arbitrarily large. A larger C" 
merely means that the distributions start to converge at 
a larger number of subsystems. Thus condition 1)12(1 re- 
quires that the energy of the total system is not concen- 
trated in only a few subsystems but distributed among 
the majority of them. Therefore, for subsystems with an 
infinite spectrum, the ratio of states where our theorem 
does not apply is negligible [l5^ . 

For applications in physics, where n is very large but 
finite, we can use the substitution 



E-Ea 
Aa 



(14) 



and write equation H13() as an integral over energies: 

Fa{E^e [Ei,E2]) = 




for all — 00 < i?i < i?2 < 00 in the limit of large n. 

The equality sign should be understood as an "asymp- 
totic" one here. The same applies to all further equations 
derived from H15|) . 

Let us finally introduce the following abbreviation for 
the density associated with the measure Pq: 



^■^i^) = /TT- A CXp I ^^2" 



27rAa 



(16) 



As a consequence the expectation value of an operator 
O, which is a function of H and thus diagonal in the 
eigenbasis, can be written 



{a\0\a) 



Wa{E)0{E)dE 



(17) 



where 0{E) is the eigenvalue of O belonging to the en- 
ergy E. If O is not a function of H , degenerate eigen- 
values of iJ, where O takes on different values, are prob- 
lematic. 

Equation H15|) can be used to estimate various quanti- 
ties of interest in physics. Among those are spectral den- 
sities and partition sums. We consider these two quanti- 
ties in the following two sections. 



IV. SPECTRAL DENSITIES 

Spectra of energy levels and thus spectral densities are 
of immense interest in the theory of quantum systems. 



They play a central role, e.g. in the analysis of chaotic 
behavior of their dynamics ja, Q ■ 

For systems, where the above theorem holds, the cal- 
culation of spectral densities is straight forward. Let us 
first consider the counting function N{E), that is the 
number of energy levels below a given threshold energy 
E. It is given by the trace of the operator Q{E — H), 
where Q is the Hcaviside step function. Since the trace of 
an operator is invariant under basis transformations, we 
choose to compute it in the basis formed by the product 
states la); 



NiE)=J2{<^mE-H)\c 
{k>} 



(18) 



where the sum extends over all states \a) of the type 
(0. According to equation (|17|l . the expectation value of 
Q{E — H) in the state \a) reads 



{a\Q{E - H)\a) ^ / Wa{E')dE\ 
Je„ 



(19) 



where Eg is th energy of the ground state of the system. 
The density of states 77 is given by the derivative of the 
counting function with respect to the energy, 'q{E) = 
dN{E)/dE; 



V{E) 



E 



Wa{E) 



E' 
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(20) 



Since the convergence of the distribution is weak, i.e. 
only on intervals of nonzero length, the derivative should 
be understood according to its definition as a linear ap- 
proximation on intervals of arbitrarily small but non van- 
ishing length. 



V. PARTITION SUMS 

The thermodynamics of a physical system is com- 
pletely determined by its partition function. It is there- 
fore of fundamental relevance to know the partition func- 
tion of the system of interest. For quantum systems its 
calculation is extremely demanding since it involves the 
complete diagonalization of the Hamiltonian. As the di- 
mension of the Hilbert space grows exponentially with 
the number of subsystems, the diagonalization of the 
Hamiltonian quickly becomes impossible even with su- 
per computers. 

Equation H15|) allows to give an analytical expression 
for the partition function at finite temperatures, which 
can easily be evaluated numerically. The partition func- 
tion is given by the trace of the operator exp(— /3iJ) with 
the inverse temperature (3. We again express it in the 
basis {\a)}: 



E(4 



-0H\ 



(21) 



The expectation values of exp{—f3H) can be computed 
using equation ((T7|l [la . llTJ I , they read 



(ale-'^^la) 



exp 



~(3Ea 



(3'Al 



erfc 



Eg - Eg + fjAl 

\/2Aa 



erfc 



V2Aa 



(22) 



where erfc(a;) is the conjugate Gaussian error function 
[l8l |. Eg is the energy of the ground state and £'„ the 
upper limit of the energy spectrum. 

Expression H22|l can be simplified further. The under- 
lying central limit theorem is valid in the limit of a very 
large number of subsystems. In that limit, the argument 
of the second conjugate error function is always much 
bigger than the argument of the first. Furthermore, it is 
always positive, which makes the second error function 
term negligible compared to the first |18j . 

Therefore, the partition function Z can be taken to 
read 



Z = 



E' 



-l3{Ea-E,) 



-erfc 
2 



Eg-Ea+ 13 Al 
V2Aa 



(23) 

where we have rescaled the energy in the first exponent, 
so that all appearing energies are positive and therefore 
exp (— /? (^Ea — Eg)) < 1. The ground state energy Eg in 
the error function is not a consequence of the rcscaling 
but stems from a cutoff in the integral (|17|l similar to the 
one in H19|) . 

In contrast to the expression for the density of states 
(|20|l . there appears one quantity in equation H23|l that 
cannot be obtained without diagonalizing the Hamilto- 
nian. This is the ground state energy Eg. 

The exact value of Eg however is not needed. The 
cutoff in H17|) at Eg is only introduced because we are 
dealing with a finite number of subsystems. In the limit 
of infinite number of subsystems, it is irrelevant since the 
Gaussian function Wa{E) decays strong enough, so that 
it becomes negligible &i E ~ Eg. A sufficiently good 
estimate of Eg can be obtained from the spectral density 

(EOll- 

However, since equation H15|l is only an approximation 
to one term in the sum (|23|l . errors due to the finite 
number of subsystems may add up. Nonetheless, the 
partition function divided by the number of states may 
be calculated instead (see section IVIII for details). 

Equations (|20|l and H23|) are the main result of the 
present paper. We now turn to verify their validity for a 
model that can be treated exactly. 



VI. NUMERICAL VERIFICATION 

In this section we present numerical tests of the two 
main results (|20|) and (|23|l for an Ising spin chain in a 



transverse field. The Hamiltonian of the chain reads 

Cn n \ 

-E<-^T<®<+1 • (24) 
i=l i=l / 

Here, af and a^ arc the Pauli matrices, 2B is the dif- 
ference between local energy levels and KB the coupling 
strength. 

The model H24|) can be diagonalized via successive 
Jordan- Wigner, Fourier and Bogoliubov transformations 
|20l |2]| . The eigenvalues of the Hamiltonian (|24|l read 



"/2 . ^. 

E^= Y^ uji ini{(p) - -j 



(25) 



l=-{n/2) + l 



where the ni{ip) are fermionic occupation numbers that 
can take on the two values and 1. The eigenfrequencies 
LUi arc given by 



UJl 



2B, K^ + 1-2K 



2ttI 



(26) 



We chose units where Planck's and Boltzmann's constant 
are equal to one, h^ ks ~ 1- 

For a finite number of spins n, a "density of states" 
can be defined with respect to certain energy bins: We 
chose the size of the bins to be B, so that the density of 
states rin{E) is defined as 

_ number of eigenstates with E^p E [E, E + B) 
VniE) = — . 

(27) 
The second quantity of interest, the partition function, 
is given by the standard expression |20| 



n/2 

Z„ = Y[ cosh (/3^) , 

/=-(n/2) + l 



(28) 



where P = T~^ is the inverse temperature. Note that 
the exact quantities, rjn{E) and Z„, carry an index n, 
reflecting the finite number of spins, in contrast to the 
values of the asymptotic approximation, 'q{E) and Z . 

Before we proceed to calculate the density of states and 
the partition function for the model H24|l with the help of 
equations (|20|1 and H23|) . let us test whether the central 
limit theorem (|15|l is applicable at all, that is whether 
conditions Hll|) and (|12|l arc satisfied. 

The energy of each spin is at least —B and at most B 
so that condition H12|) is fulfilled. The squared width A^ 
reads 



nB^ K^, 



(29) 



where n is the number of spins, and condition (|11|) is 
also met. For a large number of spins, the density of 
states and the partition function of the system at hand 
can thus, indeed, be calculated via equations (|20|) and 
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FIG. 1; Density of states for a chain of 10 spins with B = 
K — 1. The dots show the exact density rjn and the line the 
approximation 77 {rjn and 77 are defined in equations 1)27^ and 
H2()|l respectively). 
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FIG. 2: Density of states for a chain of 15 spins with B = 
K = 1. The dots show the exact density rjn and the line the 
approximation 77 (?7„ and 77 are defined in equations 1271 and 
I2OII respectively). 



In the model at hand, Ea = Ea and A^ = const. 
Therefore, the sum over all states \a) can be trans- 
formed into a sum over all energies Ea ~ k2 B — n B, 
k = 0,1, . . . ,n, 



E 

{k>} 



E 



(30) 



Figurenshows the density of states rjn{E) and its ap- 
proximation ri{E) for a chain of 10 spins, while figure [21 
shows the same plot for a chain of 15 spins. The ap- 
proximation works well despite the still small number of 
spins; furthermore, the tendency that the approximation 
improves with increasing number of spins is evident. 
Figure O shows the partition function for a chain of 100 
spins divided by the number of states, 2^°*^, as a func- 
tion of temperature. The difference between the exact 
function Zn and the approximation Z is not visible. To 
see whether the convergence improves with the number 
of spins, we have considered the maximal difference be- 
tween Zn and Z for all temperatures. 



5{n) — max |Z„ — Z\ 



(31) 



We have found the following values: (5(10) ^ 10 ^, 




ln(Z)/n , ln(Z„)/n 



10^ 10' 
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FIG. 3: Partition function for a chain of 100 spins with B = 
K = 1 divided by the number of states, 2^"". The difference 
between 2~" Zn and 2~" Z is not visible {Zn and Z are defined 
in equations H28|l and ^2'A\ respectively). 



.5(100) - 10"® and (5(1000) - 10 
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VII. DISCUSSION AND LIMITATIONS 

The approximations of the quantities rjn{E) and Z„ 
with equations (|20|l and (|23|l face some problems that 
cannot be avoided. 

Firstly, the convergence in equation (|15|l is only weak, 
i.e. the Ihs converges to the rhs for all intervals [Ei,E2] 
of nonzero length. The convergence is not pointwise. In 
the present case this means that, for example, the trace 
of a projector on a single eigenstate P = \'p){'p\ cannot 
be approximated. The problem becomes apparent, if one 
tries to calculate the partition function for zero tempera- 
ture via (jSSJ, where only the ground state is occupied. In 
the present model, for example, this state is energetically 
separated from the other states that form a quasi contin- 
uous band. As a consequence quantum phase transitions 
|2(TI | occuring at zero temperature can not be treated with 
our approach. 

The second drawback of our approach is the follow- 
ing: Each term (alOja) for some operator O that is a 
function of H is well approximated and the accuracy in- 
creases with the number of subsystems n. On the other 
hand, the number of terms in the sums H20|) and 123|) 
increases exponentially with the number of subsystems, 
for example with 2" for the spin chain. Therefore, the 
quantities ri{E) and Z can only be in good accordance 
with r]n{E) and Z„, if both are divided by the dimension 
of the Hilbert space, i.e. the number of states \a). 

This will not always be problematic, since the errors 
each term (a|0|a) carries need not all be of the same sign 
and may thus cancel each other, as in the calculation of 
the spectral densities. 

In the calculation of the partition sum, however, there 
is the following problem: Every stable system has a finite 
minimal energy, the energy of the ground state. Nonethe- 
less the probability density p6(l is, albeit very small, 
nonzero for all energies. Therefore, one needs to intro- 
duce the cutoffs in the integral ((T7|l . As can be seen from 
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FIG. 4: Logarithm of the partition function divided by the 
number of spins for a chain of 1000 spins with B = K = 1. 
The dashed line shows the exact expression In(Z„) and the 
solid line the approximation \n{Z) {Zn and Z are defined in 
equations 12811 and 12311 respectively). 



expression ()23|l . the upper limit of the integral does not 
matter. The lower limit on the other hand matters and 
becomes increasingly relevant at low temperatures. No 
matter what lower limit of the integral we take, the error 
of the approximation H23|) always has the same sign. 

Figure 0] shows the logarithm of the exact partition 
function of the spin chain and its approximation, each 
divided by the number of spins, for a chain of 1000 spins. 
The plot is done with the exact ground state energy. For 
low temperatures (T < B), the value of the approxima- 
tion is too large and the approximation fails. Therefore, 
only the partition sum divided by the number of states 
can be accurately predicted. 



VIII. SUMMARY 

We have given analytic expressions for spectral den- 
sities and partition sums of chains of quantum systems 
with nearest neighbor interaction, which are valid in the 
limit of infinitely many subsystems. We have numerically 
evaluated these expressions for a spin chain with a finite 
number of spins and compared the results with the val- 
ues obtained by exact diagonalization. The results show 
increasing accordance with growing number of spins. 

Furthermore, we have discussed the limits of the valid- 
ity of our approach and some problems that can occur, 
as well as their effects. 

The results of this paper should provide useful tools for 
the calculation of spectral and thermodynamical quan- 
tities in many systems which are intensively studied in 
present day condensed matter physics. Among those are 
spin chains [ij and strongly correlated electrons |23| . 
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